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INTRODUCTION 

This  project  originated  in  a search  to  find  mechanisms  to  explain 
the  data  obtained  at  AFCRL^  on  the  efficacy  of  lampblack  in  reducing  the 
infrared  (2-5  pm)  signature  of  aircraft  jet  engines.  These  data  were  taken 
primarily  on  stationary  aircraft  with  a variety  of  lampblack  injection  rates. 
While  the  tail-on  results  were  readily  explained  as  due  to  simple  absorption 
by  the  lampblack,  the  range  of  injection  rates  and  aspect  angles  considered 
appeared  to  require  unusually  high  absorption  coefficients  for  carbon.  Be- 
cause the  cloaking  effect  was  effective  at  all  angles,  it  did  not  seem  likely 
that  the  effect  was  due  to  scatter,  since  radiation  scattered  from  one  direc- 
tion would  then  appear  at  other  directions,  so  there  would  be  "hot"  regions. 
Cooling  because  of  the  sheer  mass  of  carbon  was  not  a sensible  solution 
because  of  the  variety  of  mass  injection  rates  over  which  the  effect  operated 
and  because  of  the  disparity  in  the  carbon-to-fuel  mass  injection  ratios 
encountered  in  military  aircraft. 

One  promising  candidate  was  a radiationless  de-excitation  (quenching) 
of  CO2  (and  also  ^0)  molecules  on  the  surface  of  the  carbon  particles.  It 
was  assumed  that  the  lampblack  mixed  with  the  hot  outer  gases  of  the  exhaust 
plume  in  a thin  shrouding  layer  surrounding  the  hottest  portion  of  the  exhaust. 
In  this  region,  the  CO2  molecules  would  interact  by  collision  with  the  surface 
of  the  carbon  particles,  losing  their  excitation  energy  to  internal  and  motional 
degrees  of  freedom  of  the  particles.  If  this  de-excitation  took  place  rapidly, 
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many  CC^  molecules  would  be  de-excited  on  each  particle.  The  particle  flow 
would  carry  the  energy  away  rapidly  to  be  dissipated  in  the  atmosphere.  At 
the  same  time,  the  de-excited  CC^  molecules  would  behave  like  a resonant 
absorber,  removing  radiative  energy  and  delivering  it  to  the  carbon.  The 
net  effect  of  this  cyclic  process  would  presumably  be  an  apparently  enhanced 
absorption  coefficient  for  the  carbon  particles. 

The  objects  of  the  project  are 

• to  examine  the  quenching  process  from  a theoretical 
viewpoint  to  obtain  estimates  of  de-excitation  rates 
to  be  expected. 

• to  conduct  a simple  experiment  to  measure  the  extinction 
coefficient  of  lampblacks  in  a static  matrix. 

• to  calculate  the  expected  signature  of  a real  jet  aircraft, 

in  flight,  with  and  without  the  quenching  mechanism  activated. 
From  the  results  of  these  three  tasks  a conclusion  is  to  be  reached  on  the 
significance  of  the  mechanism  to  infrared  signature  suppression. 

In  the  following  section  of  the  report,  the  theoretical  analysis  of 
the  quenching  of  CC^  and  1^0  molecules  on  carbon  surfaces  is  presented.  Next 
the  experimental  characterization  of  several  typical  lampblacks  is  described. 
Then  the  signature  calculations  and  an  analytical  model  are  presented.  Con- 
cluding comments  on  the  mechanism  and  its  significance  are  given  in  the  final 
section. 
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QUENCHING  OF  INFRA-REF  EMISSIONS  OF  CO,,, 
h2o  BY  COLLISIONS  WITH  GRAPHITE 


Summary 


We  have  performed  computer  simulations  for  the  collisions  of 
vibrationally  excited  CC^  and  ^0  with  a graphite  surface.  A model  poten- 
tial was  developed  to  describe  the  interaction  between  C02>  H2°  anc*  t*ie 
basal  plane  of  graphite.  This  potential  was  calibrated  against  known  bond 
strengths  and  other  properties  for  the  isolated  molecules  (CO^ , ^0)  and 
fragements  (CO,  02 , H2 , OH).  Data  for  the  interaction  of  the  molecules, 
fragments,  and  individual  atoms  (C,  H,  0)  interacting  with  the  basal  plane 
of  graphite  were  also  utilized.  In  so  far  as  these  bonding  properties  are 
accurately  known,  our  model  potential  reproduces  them. 

Using  this  model  potential,  classical  trajectories  were  calculated 
to  simulate  collisions  of  vibrationally  excited  CO 2 and  H20  with  the  basal 
plane  of  graphite.  The  basal  plane  is  not  usually  the  most  reactive  chem- 
ical site  on  a graphite  particle,  but,  statistically,  is  usually  the  most 
abundantly  exposed  plane  and,  hence,  will  collide  most  frequently  with  gas 
molecules.  However,  the  edge  planes  seem  to  bind  adsorbed  particles  more 
tightly  than  the  basal  plane.  Restricting  attention  to  the  basal  plane 
would  tend  to  underestimate  quenching. 

Since  the  interaction  potential  is  quite  realistic,  the  trajec- 
tories are  rather  long  to  compute.  As  a result,  only  a handful  could  be 
studied  for  the  various  gas  conditions,  surface  temperature,  and  other 
collision  parameters.  However,  the  results  uniformly  predict  that  when  a 
vibrationally  excited  CO 2 or  H20  (i.e.,  a potential  IR  emitter)  collides 


with  a graphite  surface,  it  will  be  destroyed  as  an  emitter.  While  there 
is  some  variation  in  the  actual  outcome  of  the  collisions,  the  one  result 
not  observed  for  any  of  the  cases  studied  is  elastic  rebound  (with  the  possi- 
bility of  subsequent  IR  emission).  Given  the  uncertainties  in  sample  size 
our  best  estimates  for  collisional  quenching  probability  is  0. 8-1.0  inde- 
pendent of  gas  or  surface  temperatures . 

In  a realistic  case,  the  graphite  surface  will  be  covered  to  an 
unknown  extent  with  adsorbed  contaminants.  Other  work  in  gas  solid  energy 
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transfer  has  shown  that  adsorbed  species  with  vibrational  characteristics 
similar  to  the  gas  molecules  (i.e.,  a near-resonant  energy  transfer  process)  will 
lead  to  enhanced  quenching.  Thus,  we  believe  that  the  simplifications  in  our 
model  will  tend  to  underestimate  quenching,  lending  further  credibility  to  our 
conclusions . 


Potential:  General  Discussion 

The  details  of  this  model  potential  for  the  interaction  of  tri- 

2 

atomic  molecules  with  a solid  surface  are  given  in  our  previous  work.  We 
briefly  review  the  appropriate  formulas  here.  We  are  interested  in  the 
interaction  of  H^O  of  CC>2  with  the  basal  plane  of  graphite.  While  probably 
not  the  most  reactive  face  of  graphite,  it  is  usually  exposed  to  a greater 
extant  than  any  other.  Hence,  it  is  most  likely  to  undergo  collisions 
with  gas  molecules.  As  indicated  below,  our  results  indicate  that  the 
basal  plane  is  a very  efficient  quencher.  One  would  expect  that  more 
reactive  sites  (e.g.  edge  planes)  would  be  certainly  no  less  efficient  and 
probably  more  efficient  due  to  the  larger  bond  strengths.  Following  our 
previous  work,  we  use  a London-Eyring-Polanyi-Sato  type  potential  function 
for  the  gas-solid  interaction.  This  potential  function,  for  three  atoms 
interacting  with  a solid  surface,  is  given  by  the  five  linear  independent 
valence  bond  structures  that  correspond  to  maximum  bonding  in  the  ground 
state  of  a six-electron  system  (see  Figure  la).  The  energy  of  the  system 
is  given  by  the  generalized  eigenvalue  equation 


HC  = SCE 


(1) 


where  the  matrix  elements  H_,  S_  are  between  the  bond  eigenfunctions 
_I^,.  . . shown  in  Figure  la.  There  are  just  six  two-body  interactions 
that  are  relevant  to  the  triatom-solid  surface  problem:  three  atom-atom 
interations  and  three  atom-surface  interations.  These  are  shown  in  Figure  lb. 
All  other  interactions  that  occur  in  the  six-electron  problem  are  set  to  zero 
here.  This  implies  a model  of  the  solid  in  which,  as  an  extended  source  of 
the  electrons  required  for  bonding, it  can  provide  the  necessary  electrons 
at  any  point  or  points  on  the  surface.  Assuming  that  the  one-electron 
orbitals  of  the  system  are  mutually  orthogonal,  we  obtain  the  matrix 
elements 


Hu  = 8Q  + 8[  (ab)  + (cd)]  - 4[  (be)  + (ac)  + (af)  + (be)] 
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H22  = + 8[ (be)  + (af)]  - 4[  (ab)  + (ac)  + (be)  + (cd)] 

H33  = 8Q  + 8 (ab)  - 4[ (be)  + (ac)  + (af)  + (be)  + (cd)]  (2) 

= 8Q  + 8(bc)  - 4[  (ab)  + (ac)  + (af)  + (be)  + (cd)] 

H55  = 8Q  + 8[ (af)  + (be)  + (cd)]  - 4[ (ab)  + (be)  + (ac)] 

H12  = H34  = h35  = H45  = 2Q  + 2[ (ab) + (be) + (af) + (be) + (cd)]  -4(ac) 

H13  = -4Q  - 4[ (ab)  + (cd)]  + 2[ (be)  + (ac)  + (af)  + (be)] 

Hl4  = -4Q  - 4[ (ab)  + (be)  -2(ac)  + (cd)]  + 2[  (af)  + (be)] 

H15  = -4Q  - 4[  (ab)  + (af)  + (be)  + (cd)]  + 2[ (be)  + (ac)] 

H23  = "4Q  ' 4[(ab^  + (bc>  “2  (ac)  + (af)]  + 2[(be)  + (cd)] 

H24  = -4Q  - 4Kbc)  + (af)]  + 2[  (be)  + (cd)  + (ab)  + (ac)] 

H25  = "4Q  ‘ 41 (bc)  + (af)  + (be)  + (cd)]  + 2t (ab)  + (ac)] 

where  Q is  the  Coulomb  integral 

Q = (abedef | H| abedef ) (3) 

and  (ab)  etc.  are  the  exchange  integrals.  The  overlap  matrix  elements 
are  given  by 


511  " S22  ~ S33  = S44  = S55  = 

512  = S34  = S35  " S45  = 2 


S13  " S14  ~ S15  = S23  " S24  S25 


- 4 


(4) 


It  is  supposed  that  the  Coulomb  integral  Q can  be  adequately  represented 

as  a sum  of  two-body  Coulomb  integrals  Q = S Q,  ,,  and  that  the  exchange 

i j 

integrals  may  be  written  as  two-body  exchange  integrals.  As  in  the  usual 
LEPS  treatment  we  assume  that  atom-atom  and  atom-surface  interactions  can 
be  described  by  Morse  potentials.  We  obtain 
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Qij  = 4~(1 +A~'~)  f(3  + Aij}  exP  ("2cyi j <ri j " ri j> > 
(2+6Aip  exp  (“ai j (ri j “ ri j>  > ^ 


(ij)  = 4(f+A~ 7)  °ij  {(1+3Aij>  exP  <-2«ij<rij-rJj>J 


(6  + 2A..)  exp  (-o'.  . (r . . - r.  .))} 
xj  r ij  xj  ij  J 


where  D , q’.jj  > rij’  an<*  ^ij  are  t^e  dissociation  energy,  Morse  parameter, 
equilibrium  distance,  and  Sato  parameter  for  the  i-j  interaction,  and 
r_  is  the  i-j  distance.  Thus,  the  matrix  elements  in  equation  (1)  are 
determined  and  hence  the  energy  E.  The  ground-state  energy  is  the  lowest 
eigenvalue  of  equation  (1). 

To  study  the  dynamics  of  gas -solid  interactions,  we  use  the 
classical  trajectory  method  of  stepwise  integration  of  the  classical 
equations  of  motion.  For  this,  we  require  the  derivatives  of  the  energy 
with  respect  to  the  three  coordinates  of  each  atom.  These  derivatives  are 
obtained  from  equation  (1)  by  the  matrix  Hellmann -Feynman  theorem: 


SE  , . dH 

iR“  ^ £ 


2^  is  obtained  by  differentiating  each  of  the  matrix  elements  of  H.  The 

eigenvectors  C are  obtained  in  the  solution  of  the  generalized  eigenvalue 
~ dE 

equation  (1)  and  so  the  derivatives  ^ and  hence  the  forces  on  each  atom 
are  easily  evaluated  at  each  step  in  the  trajectory. 

This  potential  assumes  a model  of  the  surface  in  which  the  atoms 
of  the  solid  are  fixed  at  their  equilibrium  positions  and  are  not  allowed 
to  move.  No  energy  transfer  between  the  gas  atoms  and  the  solid  surface 
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is  permitted  and  energy  terms  associated  with  geometry  changes  in  the 
substrate  atoms  are  neglected.  For  the  interaction  of  1^0  and  CO2  with 
graphite,  we  are  interested  in  possible  energy  transfer  to  the  solid  and 
so  we  must  relax  this  rigid  surface  restrictions.  The  detailed  generali- 
zation of  this  model  to  include  gas-solid  energy  transfer  was  effected 
under  a separate  AFOSR  contract  (F49620-77-C-0004) . The  results  have  been 
reported  elsewhere  and  will  be  briefly  summarized  here.  Using  the  same 
basic  potential,  VLEPS’  as  described  above,  we  include  additional  terms 
to  account  for  the  motion  of  the  solid  atoms.  These  correction  terms 


We  have 


V and  V 
r c 


(1)  V , the  restoring  force  on  each  atom  of  the  solid, 
tending  to  return  it  to  its  lattice  site. 

(2)  V , to  account  for  the  change  in  the  gas-surface 

corr 

interaction  when  the  solid  atom  is  displaced  from 
its  lattice  site. 


V = V + V + V ( 

total  LEPS  r corr 

must  vanish  when  ail  the  solid  atoms  occupy  their  lattice 


positions.  We  have  used  pair-potentials  for  and  V to  correct 
approximately  to  allow  for  the  motion  of  the  solid  atoms.  Assuming  that 
V binds  the  solid  atom  to  its  lattice  site  via  a harmonic  potential  con- 
stitutes an  Einstein  model  of  the  solid.  This  was  used  in  the  present  work. 

This  model  does  not  include  the  non-additive  corrections  for  the 
solid-solid  interactions,  nor  does  it  include  the  non-additive  corrections 
to  V . For  the  restoring  force  V we  have  used  a harmonic  potential 
binding  an  atom  of  the  solid  to  its  fixed  lattice  site.  To  account  for 
the  change  in  the  gas-solid  potential  due  to  displacement  of  the  surface 
atom  from  its  lattice  site,  we  have  used  a pair-potential  for  vcorr  con- 
necting each  gas  atom  with  each  of  the  surface  atoms  that  are  allowed  to 
move.  VT„_„  already  contains  the  interaction  of  each  gas  atom  with  each 
lattice  site,  so  to  avoid  including  it  twice,  this  lattice  site  inter- 


action must  be  subtracted  from  V 


We  take 
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corr 


(8) 


V 


Ng  Ns 

■ f ? (w(Vs>  • "“w 


where  N is  the  number  of  gas  atoms,  N is  the  number  of  solid  atoms  that 

O 

are  free  to  move,  Rg_s  is  the  distance  between  the  gas  atom  and  the  solid 
atom,  and  Rg_^  is  the  distance  between  the  gas  atom  and  the  lattice  site 
of  the  solid  atom.  As  required,  if  the  surface  atom  occupies  its  lattice 
site 


= R , 
g-s  g-l 


and  V = 0 
corr 


The  total  gas-solid  potential  for  the  moving  surface  case  is 


total 


LEPS 


Ns  o Ng  Ns 

+ Sv„  + S E 


0 s 


[«<YS>  • w(Rg-i)1 


(9) 


g 

where  VR  is  the  harmonic  restoring  force  for  solid  atoms  s. 

In  summary,  the  triatom-solid  surface  potential  that  we  have  used 
for  the  study  of  C02 -graphite  and  H2 O-graphite  interactions  includes 

(1)  all  forces,  pairwise  and  non-pairwise,  to  describe  the 
interaction  of  a gas  molecule  with  a surface  when  the 
surface  atoms  are  fixed  at  their  lattice  sites. 

(2)  pairwise  corrections  for  the  motion  of  the  solid  atoms 
away  from  their  lattice  sites. 

(3)  pairwise  corrections  to  account  for  the  change  in  the 
gas-solid  potential  due  to  displacement  of  the  surface 
atoms  from  their  lattice  sites. 

Potential:  C0„  + Graphite,  H20  + Graphite 

To  use  the  potential  described  above  for  H20  and  CC>2  interacting 
with  the  basal  plane  of  graphite,  we  require  Morse  and  Sato  parameters  for 
each  of  the  atom-atom  and  atom-solid  interactions.  The  values  of  the  Morse 
parameters  used  for  the  H-H,  H-0,  C-0,  an(j  0-0  interactions  are  readily 
available  from  spectroscopic  data  and  are  given  in  Table  I.  The  Sato  para- 
meters A for  C-0  and  0-0  in  on  were  chosen  to  be  0.3  and  0.1  respectively. 
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TABLE  I 

Morse  parameters  (in  atomic  units)  for  diatomic 
interactions . 


0-H 

H-H 

C-0 

0-0 

Do 

0.1697 

0.1744 

0.4132 

0.1903 

a 

1.2119 

1.0277 

1.2170 

1.4100 

ro 

1.8341 

1.4016 

2.1320 

2.2820 

10 


For  H^O  the  three  atom  LEPS  potential  was  that  supplied  by  A.  Gauss  (U.B.  Army 
Ballistic  Research  Laboratories)  which  the  Sato  parameters  are  equal  and  are 
given  by 

A = 100(sin6  ("^)  + 3 sin^  f-| j)  sin16  OtiMR^)^  e * 2 (10) 


where  R^  and  R2  are  the  0-H  distances  (in  a.u.),  9 is  the  angle  HCH  and 
t = arctan  (R^/R^.  This  A function  gives  rise  to  an  H2O  LEPS  potential 
in  which  the  equilibrium  geometry  and  binding  energy  are  close  to  the 
experimental  values. 

For  the  atom-solid  interactions,  the  Morse  parameters  are  chosen 
to  be  functions  of  x and  y,  where  the  x-y  plane  is  parallel  to  the  plane 
of  the  graphite  surface.  To  reflect  the  symmetry  of  the  basal  plane  of 
graphite,  we  chose  the  Morse  parameters  at  a point  (x,y,z)  to  be 


Dq  = a + 3 [cos  (3Rgl)  + cos  (9Ra2)  + cos  (3Ra3)  + cos  (BRfal)  + 
cos  (3Rfe2)  + cos  (3Rb3)  + cos  (gR^cos  (3Rbl)  + 

cos  (gRa2)cos  (3Rb2)  + cos  (gRa3)cos  (3Rb3)}  +3  {[cos2  (gRfll)-l] 
[cos2  ORbl)-l]  + [cos2  ORa2)-l][cos2  ORb2) -1]  + [cos2  (gR^-l] 
[cos2  (gRb3)-l]}  , 

and 


r = d + -§  {cos  (gR  ,)  + cos  (gR  ) + cos  (gR  _)  + cos  (gR,  ,)  + 
o j ai  az  ai  oi 

cos  (3Rb2)  + cos  (3Rb3)  + cos  (gRal)cos  (gR^)  + (12) 

cos  (gRa2 ) cos  ORb2)  + cos  (gRa3)cos  (gRb3)}  +3  {[cos2  (gRal)-l] 
[cos2  ©Rbl)-1]  + [cos2  (gRa2)-l][cos2  (gRb2)-l]  + [cos2  (BRa3)-ll 
[cos2  (0  Rb  3 ) - 1 ] } 

where 


X 

R - — 

R — — — — 4-  v 

a 1 = _Y 

■ Jz 

a3  " V3 

b2  ‘ V3  y 

X 

a2  = V3 

- y 

R - — 
bl  ~ ^3 

_ X 

^3  • ? • 

(13) 

11 


and 


0 


2n 

SV3 


where  s is  the  graphite  nearest  neighbor  bond  distance  (=  2.68  a.u.). 

The  parameters  a,  b,  c,  d,  e,  and  f are  chosen  for  each  atom-graphite  inter- 
action by  fitting  the  above  functional  forms  to  binding  energies  and  equili- 
brium heights  above  the  surface  for  each  atom  in  bridge,  on-top,  and  hole 
sites  above  the  surface.  The  values  used  were  taken  from  CNDO  calculations 
of  Messmer^and  of  Hayns ^ with  the  binding  energies  scaled  to  0.2  of  the  cal- 
culated values  as  they  suggest.  Figure  2 shows  the  different  sites  of 
approach  used;  Table  II  gives  the  scaled  values  of  Messmer  and  Haynes  as 
used  in  these  calculations.  The  resulting  parameters  are  given  in  Table  III. 
The  Morse  parameter  a for  the  atom-graphite  interactions  is  determined  by 
fitting  the  potential  energy  curve  at  one  of  the  binding  sites  to  a para- 
bola in  the  z direction;  this  gives 


O' 


= r a° 

v DQ(a.u.) 


(14) 


where  Q'q=  0.10585  for  H,  0.12603  for  0,  and  0.1392  for  C on  graphite. 

The  Sato  parameters  for  the  atom-graphite  interactions  were 
chosen  to  be  0.0  for  H and  0 in  ^0,  -0.1  for  0 in  C02,  and  -0.1  for  C 
in  C02  on  graphite.  Figures  3-5  are  equipotential  contour  plots  of 
C02  and  H20  approaching  the  graphite  surface  with  the  parameter  values 
as  given  above  and  with  the  surface  atoms  in  their  equilibrium  positions. 

The  restoring  force  for  the  solid  atoms  was  chosen  to  be  a har- 
monic oscillator  with  force  constant  determined  by  the  temperature  of  the 
surface.  The  gas-solid  atom  and  gas-lattice  site  potentials  in  the  cor- 
rection term  V rr  were  chosen  as  exponential  repulsion  terms 


W = D e 


-o'(r-ro) 


(15) 


where  D,  <y,  and  ro  were  chosen  to  be  the  same  as  in  the  gas-solid  Morse 
potential  for  the  atom  above  an  on-top  site  of  the  solid. 


Classical  Trajectories:  Methods 

With  an  approximate  potential  surface,  we  can  now  integrate 
classical  equations  of  motion  for  the  collision  of  the  triatomic  molecules 
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TABLE  II 


Binding  energies  (in  eV)  and  equilibrium  heights 
above  the  surface  (in  X)  at  different  sites  on 


the  graphite  surface. 

Bridge 


Energy  0.85 

Distance  1.0 


Energy  2.59 

Distance  1.2 


Energy  4.08 

Distance  1.2 


On-Top Hole 

H - Graphite 

0.82  0.11 

1.2  0.5 

0 - Graphite 

1.82  1.09 

1.5  1.0 

C - Graphite 

4.52  6.53 

1.3  0.6 


TABLE  III 


Parameter  values  (in  a.u.)  for  atom-graphite  interactions. 


H --  Graphite  0 - Graphite  C - Graphite 


a 

0.0244 

0.0814 

0.1846 

b 

-0.0068 

-0.0138 

0.0185 

c 

0.0011 

-0.0442 

-0.0370 

d 

1.6535 

2.2393 

2.0790 

e 

-0.2362 

-0.1228 

-0.2830 

f 

0.7770 

0.7937 

0.4610 

Jl 
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Figure  3 a.  Equipotential  contours  (in  eV)  for  an  atom  approaching 
surface  of  graphite  in  (X,Z)  plane  as  defined  in 
Figure  2 (y =0.0)  (a)  H atom. 


Figure  5a. 


Equipotential  contours  (in  eV)  for  H2  approaching  the 
graphite  surface  in  Figure  2.  The  H2  bond  is  held 
parallel  to  the  surface  and  the  bond  midpoint  is 
fixed  over  site  "A"  (7A). 


-*  ZZ.  'WZ f&*~ 


with  the  surface.  Observing  the  outcome  of  such  trajectories  will  provide 
an  estimate  of  the  quenching  efficiency.  The  initial  conditions  for  the 
AI^  - graphite  trajectories  are  specified  as  follov’s.  A given  trajectory 
has  preselected  vibrational  quantum  numbers  (v^,  V2,  v^),  rotational  energy, 
translational  energy,  direction  of  approach  of  the  molecule  to  the  surface 
and  temperature  of  the  surface.  The  initial  geometry  of  the  molecule  is 
selected  from  the  inner  and  outer  turning  points  for  the  AB  bond  distances 
and  the  angle  BAB  determined  from  a normal  mode  analysis.  This  gives  a total 
of  eight  possible  geometries  for  the  molecule  and  each  geometry  is  selected 
in  turn  by  successive  trajectories.  A variable  vibrational  phase  is  intro- 
duced by  setting  the  initial  z coordinate  of  the  center  of  mass  of  the 
molecule  to  be 


Z = A + T.V  . X (16) 

c.m.  z 

where  A is  a constant  (taken  to  be  10.0  a.u.  for  CO2  + graphite  trajec- 
tories and  8.0  a.u.  for  1^0  + graphite  because  of  the  smaller  step  size 
used  and  therefore  longer  integration  time.  V2  is  the  z component  of 
velocity  of  the  center  of  mass  of  the  molecule  and  x is  a pseudorandom 
number  between  0 and  1.  T is  a classical  half-period  of  vibration  of  the 
molecule  - that  of  the  smallest  vibrational  frequency  (1033.8  a.u.  for 
CO2  and  432.3  a.u.  for  1^0) . The  velocity  of  the  center  of  mass  is 
determined  by  the  chosen  speed  and  the  direction  of  approach  of  the  mole- 
cule given  by  angles  9,  0 in  the  usual  notation.  0 was  chosen  to  be 
0.0  for  these  calculations.  The  initial  position  of  the  center  of  masss 
is  determined  by  Zcm,  0,  0 and  an  aiming  point  on  the  graphite  surface 
which  is  the  point  on  the  surface  that  the  center  of  mass  of  the  molecule 
would  pass  through  if  no  interaction  with  the  surface  were  to  take  place. 
For  the  trajectories  computed  here,  the  aiming  point  is  chosen  to  be  one 
of  the  four  sites  shown  in  Figure  2. 

The  initial  orientation  of  the  molecule  is  specified  by  the 
usual  Euler  angles  9',  0',  Y between  0 and  2n  and  cos  9'  between  -1 
and  1.  The  intial  velocity  of  the  individual  atoms  is  then  determined 
from  the  velocity  of  the  center  of  mass  of  the  molecule,  since  zero 
rotational  energy  was  chosen  for  these  calculations. 


For  the  graphite  surface  atoms  that  are  allowed  to  move,  the 
initial  positions  are  chosen  to  be  the  equilibrium  positions  of  the  atoms 
(at  the  lattice  sites).  The  Einstein  temperature  of  the  graphite  surface 
is 


and 


surface 

0 

Einstein 


o 0 • 6 9 


bulk 

Debye 


(17) 


bulk 

0 


= 420°K 


Debye 

80  that  surface  0 <18> 

0 = 217  K 

Einstein 

This  leads  to  a harmonic  force  constant  for  the  graphite  solid  atoms  of 
0.0103  a.u.  The  initial  kinetic  energy  of  each  solid  atom  was  chosen 

equal  to  the  average  energy  of  a harmonic  oscillator  at  the  given  surface 

o 

temperature  and  with  an  Einstein  temperature  of  217  K.  Kinetic  energies 
of  the  solid  atoms  for  surface  temperatures  of  500°K,  1000°K,  and  2000°K 
are  given  in  Table  IV. 

The  step  sizes  used  in  the  integration  procedure  were  3.5  a.u.  for 
CO2  on  graphite  and  0.7  for  1^0  on  graphite,  the  difference  being  due  to 
the  form  of  potential  function  used.  The  functional  form  used  for  A in 
the  1^0  potential  leads  to  instability  of  the  integration  unless  the 
smaller  step  size  is  used.  To  end  a trajectory,  the  usual  end  conditions 
for  a three  atoms-surface  trajectory  were  checked;  in  addition,  a tra- 
jectory was  considered  to  be  complete  if  the  number  of  integration  steps 
taken  exceeded  a given  maximum  (30,000  steps  for  CO2  on  graphite  and 
40,000  steps  for  1^0  on  graphite). 


Classical  Trajectories:  Results 

The  results  of  the  1^0  and  graphite  and  CO2  and  graphite  tra- 
jectories are  given  in  Tables  V and  VI.  All  of  the  trajectories  we  were 
able  to  compute  resulted  in  the  molecule  (or  some  fragment)  remaining 
bound  to  the  surface.  For  the  1^0  and  graphite  system,  all  trajectories 
resulted  in  the  1^0  molecule  being  stuck  on  the  graphite  surface. 
in  its  (0,0,1)  vibrational  state  was  used  and  the  molecule  was  initially 
aimed  perpendicularly  towards  an  on-top  site  of  graphite;  the  solid  atom 

O 

in  that  site  was  allowed  to  move.  A surface  temperature  of  500  K was 
used  with  ^0  translational  energies  of  2000°K  (10  trajectories)  and  1000°K 
(4  trajectories).  Four  more  trajectories  were  computed  with  a graphite 
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TABLE  IV 


TABLE  V 

1^0  + Graphite  Trajectories 


V1  v2  v3  0 Aira  Polnt  Surface  Translational  Number  Number 

temperature  energy  of  of  each 

(°K)  (°K)  trajectories  outcome 


f 


TABLE  VI 


CO 


2 + Graphite  Trajectories 


V1 

V2 

V3 

6 

Aim  Point 

Surface 

Temperature 

(°K) 

Translational 

energy 

(°K) 

Number 

of 

Trajectories 

Number 
of  each 
outcome 

1 

0 

1 

0.0 

A 

500 

2000 

4 

4b 

0 

0 

1 

0.0 

A 

2000 

2000 

4 

4b 

0 

0 

1 

0.0 

A 

500 

1000 

4 

4b 

0 

0 

1 

0.0 

A 

500 

2000 

5 

5b 

0 

0 

1 

15.0 

A 

500 

2000 

2 

lb,  lc 

30.0 

2 

2b 

45.0 

2 

2b 

60.0 

2 

2b 

0 

0 

1 

0.0 

B 

500 

2000 

4 

3b,  lc 

C 

4 

4b 

D 

4 

4b 

Outcome  b corresponds  to  an  0 atom  stuck 
molecule. 


Outcome  c corresponds  to  a gas-phase  0 at 


2i 
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possible  to  calculate  only  a small  sample  of  trajectories  for  each  set  of 


initial  conditions.  The  time  taken  was  to  large  partly  because  of  the  small 
step  size  necessary  to  retain  stability  in  the  integration  and  partly  because 
all  these  trajectories  ended  only  when  the  maximum  number  of  steps  (40,000) 
had  been  taken. 

For  CO^  and  graphite,  the  trajectories  took  somewhat  less  time  to 
compute  (about  250  secs  each),  both  because  of  the  larger  step  size  used 
and  because  the  majority  of  trajectories  did  not  take  the  maximum  number 
of  steps  but  finished  when  a CO  molecule  had  left  the  surface,  leaving  an 
adsorbed  0 atom.  Trajectories  were  computed  for  both  the  (0,0,1)  and  the 
(1,0,1)  states  of  CO^ . Surface  temperatures  of  500  K and  2000  K were  used 
and  CO2  translational  energies  of  2000  K and  1000  K.  The  aiming  point  on 
the  surface  was  varied  as  was  the  angle  of  approach  0.  In  all  cases  except 
two,  the  trajectory  resulted  in  an  adsorbed  0 atom  and  a desorbed  CO  mole- 
cule. In  the  remaining  two  cases,  the  outcome  was  a gas-phase  0 atom  and 
an  adsorbed  CO  molecule. 
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MEASUREMENT  OF  THE  EXTINCTION  COEFFICIENT  FOR  CARBON 


The  object  of  this  task  of  the  project  is  to  make  a simple  measure- 
ment of  the  extinction  coefficient  of  carbon  particles,  at  densities  appropri- 
ate to  the  shroud  concept,  in  a static  matrix  of  KBr.  The  matrix  material 
used  is  a standard  spectroscopic  matrix,  well-characterized  and  readily 
available . 

To  determine  the  densities  required,  a simple,  2-core  model  of  the 
plume,  described  in  Section  IV,  was  adopted.  In  this  model,  the  plume  is 
assumed  to  consist  first  of  an  inner  core  of  hot  gases  with  base  diameter 
equal  to  that  of  the  jet  exhaust,  and  with  apex  pointing  rearward.  Surround- 
ing this  cone  is  a second  cone  of  the  same  cone  angle,  with  apex  pointing 
forward,  and  intersecting  the  first  cone  at  the  edge  of  the  exhaust  opening. 

The  particles  are  assumed  to  be  entrained  between  these  two  cones.  The 
particle  density  for  lampblack  injected  into  this  volume  were  calculated  and 
taken  as  representative. 

Lampblack  samples  obtained  from  Cabot  in  Cambridge,  Mass.,  with 
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nominal  diameter  1500  A were  used.  The  lampblack  was  mixed  with  KBr  powder 
in  the  quantities  needed  and  portions  of  the  mixture  were  pressed  into  disks 
for  spectroscopic  analysis  in  the  2.6  to  6.0  pm  wavelength  range.  Several 
disks  (5  or  6)  were  chosen  from  those  fabricated  on  the  basis  of  total  density 
variation.  All  densities  were  required  to  fall  within  10%  of  the  mean. 

For  an  aircraft  of  velocity  1000  ft/sec.  with  exhaust  diameter  of 
26.6  inches,  the  particle  density  corresponding  to  an  injection  rate  of 
R (in  oz/sec)  is 

pp  = 3.14  x 10  ^R  (gm/cm^) 

_6  3 

The  fabricated  density  was  3.4  x 10  (gm/cm  ) , corresponding  to  R = 1.1  oz/sec. 

”*6  3 

Five  standard  disks  with  density  p =0  and  five  with  density  3.4  x 10  gm/cin 
were  measured,  in  order  to  improve  the  statistics  of  the  measurements.  Normalized 
intensity  ratios  for  each  wavelength  were  obtained  by  dividing  the  intensity 
ratio  of  the  sample  average  by  that  of  the  standard  (p^  = 0)  average.  The  average 
cross  section  of  a single  particle  was  obtained  from 

o = - ln(I/IQ) /NT 
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where  I/Iq  is  the  normalized  intensity  ratio,  T is  the  mean  sample  thickness, 

and  N is  the  particle  density.  The  values  obtained  fell  in  the  range 
-11  2 

(2.5  - 8.8)  x 10  cm  /particle,  which  are  comparable  in  magnitude  to 
6 

published  values.  The  measurements  are  summarized  in  Table  VII. 


CALCULATION  OF  SHROUD  EFFECT  ON  SIGNATURE 
Plan  of  Computations 

The  original  plan  of  this  research  was  to  carry  out  a semiquantita- 
tive  analysis  of  the  effect  of  quenching  on  the  shroud  effect  of  carbon  particles. 
Briefly,  it  was  planned  to  make  rough  estimates  of  the  mixing  layer  thickness 
using  data  available  in  the  literature.  Then,  using  Air  Force  data  on  aircraft 
velocities,  exhaust  velocities,  temperature  profiles  and  CO2  and  1^0  concentra- 
tion profiles,  the  signatures  with  and  without  the  carbon  shroud  would  be  com- 
puted using  existing  radiative  transfer  programs  modified  to  include  the 
additional  effects. 

Early  in  the  program,  it  was  pointed  out  that  the  Air  Force  had 
developed,  under  contract  with  General  Electric,  a comprehensive  set  of  programs 
for  the  calculation  of  jet  engine  signatures,  and  that  the  programs  known  as 
SCORPIO  III,  were  installed  on  the  Air  Systems  Division  computer  coupler  at 
Wright-Patterson  AFB.  After  some  discussion,  it  was  decided  to  use  these 
programs  for  the  following  reasons: 

• They  were  proven  programs. 

• They  could  be  easily  modified. 

• They  were  well  documented. 

• They  included  programs  for  determining  the  flow  field. 

• They  included  all  relevant  absorption  and  scattering 
effects. 

Further,  the  authors  of  the  programs  were  located  nearby  and  were  willing 
to  cooperate  in  implementing  the  needed  changes.  These  programs  offered 
a considerable  expansion  of  the  sophistication  of  the  signature  calculations 
at  no  apparent  increase  in  cost  to  the  project.  Finally,  as  discussed  below, 
a later  version  of  the  programs,  called  SCORPIO  IIIA,  included  calculation  of 
the  flow  field  of  a carbon  shroud  and  its  effect,  without  quenching,  on  the 
IR  signature. 

There  ensued  an  unfortunate  series  of  events  which,  ultimately, 
prevented  the  use  of  these  programs  for  the  purposes  of  the  project  and 
reduced  the  available  funding  to  an  extend  that  made  it  impossible  to 
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return  to  the  original  plan.  A brief  description  of  these  difficulties  follows, 
after  which  the  model  calculations,  developed  to  replace  the  more  realistic 
computations,  will  be  described. 

Difficulties  With  Air  Force  Codes 

After  some  delays  caused  by  the  need  for  large  amounts  of  permanent 
file  space  for  the  SCORPIO  programs,  initial  attempts  were  made  to  run  the 
codes  using  the  test  data  set  installed  with  the  source  codes.  The  data  set 
was  not  accepted  by  the  code,  and  one  source  suggested  that  SCORPIO  III  had 
never  been  successfully  run  on  the  current  operating  system.  When  installed 
the  data  set  would  not  run,  the  data  set  for  the  test  case  in  the  SCORPIO  III 
Users  Manual  was  reconstructed  manually.  This  data  set  also  did  not  run  to 
completion,  and  the  results  that  were  obtained  did  not  match  the  results  in 
the  Users  Manual.  At  this  point  it  was  discovered  that  the  installed  code  had 
been  modified  and  that  it  was  an  early  version  of  a code  known  as  SCORPIO  IIIA, 
which  incorporated  shrouding  particles.  A classified  General  Electric  report 
describing  the  SCORPIO  IIIA  modifications  was  obtained.  Using  these  modifica- 
tions together  with  some  of  the  data  changes  found  in  the  originally  installed 
test  cases,  it  was  possible  to  obtain  a data  set  that  was  completely  compatible 
with  the  installed  code.  All  of  these  steps  were  performed  laboriously  by 
step-by-step  comparison  of  coded  lines,  data  items,  and  the  various  documenta- 
tions available.  When  data  assembly  was  completed  and  run,  it  was  found  that 
the  code  still  did  not  run  to  completion  because  of  a coding  error.  The  results 
that  were  obtained  prior  to  encountering  the  coding  error,  however,  now  agreed 
approximately  with  the  test  results  of  the  Users  Manual  when  the  particle 
mass  injection  rate  was  set  to  zero.  At  this  point,  it  appeared  that  the 
codes  could  be  made  to  work,  with  help  from  General  Electric  personnel.  However, 
these  efforts  had  exhausted  80%  of  the  funds  allocated  for  the  computations. 
Further  attempts  to  use  the  code  revelaed  several  other  coding  errors  which 
were  found  and  corrected  with  help  from  G.  E.  However,  it  was  never  possible 
to  obtain  a completely  successful  run  matching  all  of  the  output  reported  in 
the  Users  Manual.  And,  a fatal  flaw  still  remains — the  radiation  from  the  hot 
engine  parts  does  not  interact  with  the  plume.  This  is  a major  contribution 
to  the  signature,  and  this  interaction  is  essential  for  operation  of  the 


quenching  mechanism.  Since  the  exact  cause  of  and  the  extent  of  the  correction 
needed  for  this  problem  were  not  known,  it  was  decided  to  halt  this  portion  of 
the  program  and  examine  the  options  available.  There  appeared  to  be  three 
viable  options,  all  requiring  some  additional  funding.  These  were 

• Abandon  the  attempt  to  use  SCORPIO  IIIA  and  return 
to  the  original  plan  to  use  BCL  codes. 

• Subcontract  the  calculations  to  General  Electric  to 
carry  out  on  the  latest  version  of  SCORPIO  IIIA, 
installed  and  routinely  run  on  their  computing  equipment. 

• Abandon  large-scale  computations  and  resort  to  model 
calculations  to  give  an  indication  of  the  relevance  of 
quenching  to  signature  suppression. 

Because  of  limited  funding,  the  third  of  these  options  was  taken. 

The  experience  with  the  SCORPIO  codes  was  disappointing,  but  not 
entirely  wasteful.  The  completed  runs  do  contain  the  temperature  profiles 
plume  flow  data,  and  exit-plane  radiation  data  necessary  to  perform  quanti- 
tative calculations. 

Description  of  the  Model 

The  actual  plume  situation  is  modeled  as  follows:  VJe  have  a mixture 
of  carbon  particles,  air,  and  both  excited  and  unexcited  C0£  molecules.  Radia- 
tion is  incident  upon  this  mixture.  The  unexcited  CC>2  can  be  excited  by  the 
radiation  and  excited  CO^  can  either  radiate  into  the  field  or  can  be  de-excitod 
without  radiation  by  collision  with  the  carbon  particles.  Each  carbon  particle 
is  treated  as  a sphere  of  radius  a,  surrounded  by  a mixture  of  air  and  CO^  in  a 
spherical  region  of  radius  b.  The  CO^  molecules  can  diffuse  to  the  surface, 
where  a fraction  y of  them  are  de-excited.  The  radius  b is  chosen  to  enclose 
that  volume  of  the  air-CO,,  mixture  determined  by  the  average  density  of  carbon 
particles  in  the  plume.  At  radius  b,  the  excited  C02  concentration  is  unchang- 
ing, i.e.,  a pure  "reflection"  boundary  condition  is  used.  The  diffusion  equa- 
tion is  to  be  solved  for  this  model  to  obtain  CO^  concentration  profiles  as  a 
function  of  time.  The  total  CC>2  concentration  is  held  constant  throughout.  C02 
alone  is  treated  as  a radiating  species  because  it  is  generally  regarded  as  the 
most  important  molecule  in  signature  determination. 
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The  use  of  a diffusion  approach  in  the  calculations  is  motivated 
by  the  need  to  determine  the  time  scale  in  which  the  quenching  mechanism 
operates.  A short  time  scale  would  indicate  that  the  quenching  is  essentially 
complete  at  the  exit  plane,  while  a longer  time  scale  would  indicate  a more 
gradual  rise  in  the  quenching  enhancement . A very  long  time  scale  would  mean 
that  the  mechanism  is  ineffective,  since  the  quenched  region  would  be  far 
behind  the  aircraft.  It  is  also  important  to  realize  that  the  time  scale  on 
which  the  carbon  particles  reach  equilibrium  with  the  air  in  the  plume  (most 
of  the  plume  gas  is  nitrogen)  is  very  short.  Thus,  we  can  take  the  carbon 
particles  to  have  the  same  temperature  and  velocity  as  the  air  in  the  slip- 
stream where  mixing  is  assumed  to  occur. 

Mathematical  Description  of  the  Model 

The  diffusion  equation  for  the  model  situation  is 

3c*  n 3 c . 2D  3c  * , 

-rr  = D — 5 + — s c + Ike  (1) 

dt  . z r dr  o T 

3r 

where 

k 

c = Concentration  of  excited  C0„ 

^ * 

c,j,  = Total  concentration  of  CO^  = c + c 
D = Diffusion  coefficient  of  CC^  in  ^ (or  Air) 

I = Intensity  of  radiation  field 

k 

k = cross-section  for  c -*■  c transition 
s = X + Ik 

° * 

X = net  decay  rate  for  c -+  c 

The  boundary  conditions  to  be  imposed  are 


9c  (b,  t)  = 
8r 


3c*(a>t)  1 j](V  *,  v 
3r  4 D U,t; 


I 

I 


- 


Here 


a ■=  radius  of  carbon  particle 

b = radius  of  sphere  of  influence  for  one  particle 
Y = fraction  of  CO*  molecules  striking  carbon  surface 
which  are  de-excited  without  radiation 
v = average  CC^  molecular  velocity. 

The  derivation  of  Eq.  ( 1) — ( 3)  follows  standard  procedures.  The  initial 
condition  for  Eq.  (1)  is  written 

c*(r,o)  = c*  (4) 

where  c is  the  concentration  of  excited  C0o  in  the  plume  at  t = 0. 

The  substitution 


transforms  Eq.  (1)  into 


3u  ^ 3 u 


with  initial  condition 


* IkCTs 

“(r.o)  = (c  o - — >r 


and  boundary  conditions 


(5) 


(6) 


(7) 


(8) 


where 


3u(a, t) 
3r 


1 7kcT 

(h  + — )u(a,t)  + ha  

SL  S 


h 


1 xi 
4 D 


(9) 


(10) 


36 


Our  problem  has  now  been  reduced  to  one  analogous  to  heat  flow  in  one 
dimension  through  a rod  under  very  general  boundary  conditions: 

1.  There  is  radiation  with  conductivity  sq  along  the 
rod's  sides  into  a medium  at  u = 0; 

2.  There  is  radiation  at  4 = b with  conductivity  1/b 
into  a medium  with  u = 0; 

3.  There  is  radiation  at  r = a with  conductivity  h + 
into  a medium  with 


1 + ha  s 


We  follow  standard  methods  to  obtain  a solution.  First,  we  put 


u(r,t)  = v(r)  + w(r,t) 


Then  we  can  obtain  the  two  differential  equations 


D — ^ - s v = 0 
3r2  ° 


^ ■ -i 


3v(a)  „ . U , , . , 

-g?-  = (h  + — )v(a)  + ha  — 


3w  _ 3 w 
— = 0 — - so  w 

Jr 


* IkcT 

w(r,0)  = r (cq ^ -)  - v(r) 

o 


3w(a,t) 

3r 


(h  + — ) w ( a , t ) 


This  technique  gives  us  a simple  time-independent  equation  for  v satisfying 
the  full  boundary  conditions  and  a time  dependent  equation  for  w satisfying 
a set  of  simpler,  homogeneous,  boundary  conditions. 

The  solutibn  of  the  v(r)  equations  are  easily  obtained  by  Laplace 
transformation  methods  and  is 

halkc„  coshp  (b-r)  sinhp  (b-r) 


(h  + — - — ) coshp  Si  - [— — 
a b o bp 


PQ]  sinhyo*. 


where 


/So 


H = b-a 


The  w(r,t)  equations  are  straightforward  to  solve  by  variable  separation. 


w(r,t)  = R(r)T(t) 


Then,  we  find 


T(t)  + T e 
o 


R(r)  = Acospr  + Bsinpr 


where  a is  a separation  constant  and 


/“ 

p = v — 


The  boundary  conditions,  Eq.  (17)  and  (18),  lead  to  a consistency 
condition: 


tanpH  = 


. . 1 _ I 

w(  a b) 

2 j.  1 . In 

w + b (h  + 


which  has  real,  discrete  eigenvalues  p^.  Using  either  of  the  boundary 
conditions  allows  R(r)  to  be  written  as 


(26) 


where 


R (r)  <*  X (r) 
n n 


Xn(r)  = Nn[cosyn(b-r)  - - — sinun(b-r) ] 


and  is  a normalization  factor.  Thus,  the  general  solution  can  be 
written 


w(r,t)  = £ AX  (r)e  ^So  + n*^ 

n=l  n n 


Only  positive  eigenvalues  are  used  because  of  symmetry.  There  is  no 
eigenvalue  for  n = 0 as  long  as  h i*  0.  The  are  found  from  the  initial 
condition: 

* IkcT 

w(r,o)  = r(cQ  - )-v(r)  = £ \\(,r) 

o n=l 


The  X^  are  orthogonal,  so  we  find 

Ike 

A = Jb[(c* — — r-v(r)  ]X  (r)dr 

n ' o s n 

a o 


The  time  scale  of  the  evolution  of  c is  determined  by  the 
lowest  eigenvalue,  y^,  and  we  define  the  time  constant  of  the  system  to  be 


11  so  + D"l 


To  determine  y^,  we  have  to  solve  the  eigenvalue  equation,  Eq.  (25) 
Define 

K - V* 

a i 
£ " b < 1 

p = 1 + ha  > 1 


i. 


Then  Eq.  (25)  can  be  written  as 


tans  - fcl)<,-Ci;  2 
z.  e + (1-e)  p 


(33) 


Generally,  we  expect  b>>a,  so  e<<l.  Then  we  can  rewrite  (33)  as 


cot£  = — £ + — 

P C 


(34) 


and  we  have  a one-parameter  set  of  eigenvalues  to  examine.  Since 
£ 

— is  small,  the  solutions  will  lie  near  the  solution  of 


cote  = 1/e  or  tans  = £ 


(35) 


These  solutions  are  tabulated.  We  find  ^ 4.5  so 


„ 4.5 

U.  'Xi  

Mi  a 


(36) 


The  time  constant  can  now  be  written  as 


T1  = 


s SL2  + (4.5  D)2 
o 


(37) 


2 2 

For  small  s X,  , using  D 'v  0.5  cm  /sec  for  typical  plume  temperature  and 

° -4  -8 

pressure,  we  estimate  t..  = 10  sec  for  l = 1 mm  and  x1  = 10  sec  for 

2 1 -1  1 

l = 1 pm.  If  s H is  large,  then  t = s , independent  of  the  value  of  i. 

° 2 

We  will  see  below  that  s S,  is  small  for  our  system. 


Signature  Calculations 

Having  obtained  solutions  to  the  diffusion  equation,  we  are  in 
a position  to  set  up  the  formulae  for  the  calculation  of  IR  signatures.  It 
is  not  intended  that  such  calculations  be  complete,  but  only  that  some 
useful  estimate  of  the  enhancement  of  absorption  by  quenching  be  obtained. 
The  inputs  from  the  diffusion  equation  solutions  are  the  concentrations, 
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c and  c . The  first  step  in  a signature  calculation  is  to  calculate  the 
_ 

concentrations,  c and  c , averaged  over  the  volume  b.  We  define  this 
average  by 

c*(t)  = Ar  Jbc*(r,t)r2dr  (38) 

bJ  a 

i.e.,  the  volume  average  over  the  characteristic  volume.  This  average, 
combined  with  the  exhaust  flow  rate  and  the  air  speed  of  the  aircraft, 
determines  a variation  of  concentration,  from  exit  plane  backwards  into 
the  distant  plume,  arising  from  the  time  dependence  of  c . That  is,  as 
we  move  back  from  the  exit  plane,  we  pass  to  concentration  corresponding 
to  later  time  epochs  in  the  solution  of  the  diffusion  equation. 

The  next  step  would  normally  be  to  set  up  some  model  for  the  plume 
structure,  and  carry  out  numerical  calculations  of  the  signature  for  various 
aspect  angles.  Alternatively,  we  can  simply  calculate  the  total  extinction 
coefficient  for  a representative  situation.  This  enables  us  to  reduce  our 
assumptions  about  the  plume  to  the  bare  minimum  and  to  carry  out  a semi- 
analytic  analysis.  We  assume  that  the  plume  consists  basically  of  two 
parts,  a central  hot  core  and  a flaring-out  region  of  cooler  gases  (see 
Fig.  6).  The  cone  angle  of  both  regions  are  taken  to  be  the  same.  Then, 


Fig.  6.  Schematic  of  plume-shape  assumption. 


we  assume  the  particles  are  mixed  into  the  cool-gas  region.  The  median 
thickness  of  this  region  and  the  mean  radius  of  the  median  annular  cross- 
section  are  both  equal  to  d^ . If  v is  the  aircraft  velocity,  then  new 
volume  is  swept  out  at  the  rate 


* 1 .2 

V = tt  dj  v 


(39) 


If  particles  are  injected  at  the  rate  of  R oz/sec,  then  the  particle  mass 
density  in  the  shroud  is 


p = R/V  = (£)3  p 
p be 


(40) 


since  the  radius  of  a particle  is  a and  each  particle  has  an  associated 

volume  corresponding  to  a sphere  of  radius  b.  In  (40),  p is  the  density 

3 C 

of  graphite,  p^  = 2.25  gm/cm  . Using  the  data  from  Table  VIII,  we  obtain 
from  Eq.  (39)  and  (40) 

b/a  = 104/R1/3  (41) 


In  order  to  cover  the  practical  range  of  commercially  available 
lampblacks  and  injection  rates  that  are  reasonable,  we  consider  the 
values 

a = 500,  1000,  5000,  10,000  A 
R = 1,  10,  100  oz/sec. 

To  determine  the  time  evolution  of  the  quenching  effect,  we  use 

the  values  from  Table  VIII  to  get  IR  = 0.213/sec.  Thus,  since  X = 500/sec, 

2 2 

sq  is  dominated  by  X.  However,  sq  is  still  small  compared  to  (4.5  D)  in 
Eq.  (37) . Values  of  are  given  (in  microseconds)  in  Table  IX.  For 
v * 1000  ft/sec,  the  quenching  effect  is  essentially  complete  for  distances 
larger  than  about  1 cm  behind  the  exit  plane.  Hence,  diffusion  does  not 
appear  to  create  any  significant  slowing  of  the  quenching  mechanism. 

The  next  step  is  to  examine  the  effective  total  extinction  coeffi- 
cient. Light  passing  through  the  shroud  is  attenuated  both  by  absorption  on 
the  lampblack  particles  and  by  absorption  through  excitation  of  the  un-excited 
(X>2  molecules.  The  carbon  contributes  a term 


From  SCORPIO  III  Codes: 

dj  = 26.6  in  = 67.6  cm 

A 

V = 1000  ft/sec.  = 3.1  x 10  cm/sec 


P = 0.45  Atm. 

T = 1230°R  = 328°K 


at  x 

o 


Radiant  Intensity  = 40  Watts/steradian 

-3  2 17 

I = 7 x 10  Watts/cm  = 1.5  x 10  photons/sec 
o 

(at  x , Y and  4.3  pm  wavelength  band) 
o o 
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c^,  = 5 x 10  molecules  of  CO^/cm 


From  literature  and  direct  calculation: 

X = 500  sec  ^ 

—18  2 

k = 1.4  x 10  cm  /molecule 

2 

D for  CO2  in  Nj  = 0.165  cm  /sec 

2 

AT  P,T  above,  D = 0.43  cm  /sec 

3 

Pc  = density  fo  Carbon  = 2.25  gm/cm 

v = thermal  velocity  for  CO^  molecule  at  328°K 
4 

= 3.81  x 10  cm/ sec. 


ac  = ac  x concentration  of  carbon  particles 


= 0 /4  tt  b3)  = 0.239  — 
c 3 b3 


to  the  extinction  coefficient.  Using  the  largest  measured  value  from  the 
third  section  of  this  report,  we  find 


-11  -3  -1 

a = 2.1  x 10  b (cm  ) 
c 


Values  of  a are  listed  in  Table  IX. 
c 

The  contribution  from  unexcited  (X^  is  obtained  by  calculating 
the  long-time,  mean  concentration  of  unexcited  CO2  in  the  plume.  To  do 
this,  we  need  to  know  the  intensity  throughout  the  thickness  of  the  shroud. 
Examination  of  the  long-time  solution,  from  Eq.  (5),  (11),  and  (28) 


c (r,«)  = 


+ v(r)/r 


with  v(r)  given  in  Eq.  (19),  shows  that  c*  is  a nonlinear  function  of  I. 
However,  since  A is  the  dominant  factor  in  so,  this  nonlinearity  is 
effectively  eliminated. 

The  long-time  behavior  of  c*  is  strongly  dependent  on  the 
parameter  pQ,  given  in  Eq.  (2).  For  the  parameter  values  of  significance 
to  us,  we  have 


/ o ...  -1 

PQ  = / — = 34.1  cm 


The  arguments  of  the  hyperbolic  functions  are  at  most  pQb  in  magnitude. 
Referring  to  Table  VIII,  we  see  that  the  largest  value  of  p^b  is  0.355, 
occurring  for  a = 1 pm,  R = 1 oz/sec.  Considering  the  smallness  of  the 
pQb,  expansion  of  the  expression  for  c*  in  powers  of  pQb  appears  to  be  one 
way  to  simplify  the  analysis.  The  interrelationship  among  the  parameters 
in  the  full  expression,  is  very  delicate,  so  we  must  proceed  cautiously, 
c*  is  given  by 


Table  IX.  System  Parameters  and  Calculated  Results 


a 

pm 

b 

pm 

T1 

p sec 

^o15 

a 

c 

(cm“l) 

S/ 

D 

P 

max 

R = 

1 

02/ sec : 

0.05 

5.20 

0.071 

0.0177 

0.149 

0.882 

0.414 

0.10 

10.4 

0.283 

0.0356 

0.0187 

0.776 

2.91 

0.50 

52.0 

7.06 

0.177 

1.49  x 

10-4 

0.303 

142.1 

1.00 

104.0 

27.9 

0.354 

1.87  x 

10-5 

0.131 

492.4 

R = 

10 

02/ sec: 

0.05 

2.42 

0.0149 

8.21  x 

10 

-3  1.49 

0.987 

0.0464 

0.10 

4.83 

0.060 

0.0164 

0.186 

0.972 

0.365 

0.50 

24.2 

1.49 

0.0821 

1.49  x 

10-3 

0.813 

3.817 

1.00 

48.3 

5.96 

0.164 

1.86  x 

10-4 

0.598 

224.8 

R = 100  02/ 

sec : 

0.05 

1.12 

3.1  x 10-3 

3.81  x 10-3 

14.9 

0.999 

4.70  x 10-3 

0.10 

2.24 

0.0122 

7.63  x 10-3 

1.86 

0.997 

0.0376 

0.50 

11.2 

0.307 

0.0381 

0.0149 

0.977 

4.60 

1.00 

22.4 

1.23 

0.0763 

1.86  x 10-3 

0.937 

35.3 

c*(r.o)-c*(r1h)  coshyo(b-r)-  — g sinhuo(b-r) 

7T — = "(ha)  2 z , 

c (r1o)  r (1+ha-  t-)  coshy  l r(l+ha-  -r-(y  b)  )sinhy  l 


° hDb 


where 

* IkCT 

c (r,o)  = — 

o 

and  h Is  given  by  Eq.  (10).  Using  = 0.8  and  v,  D from  Table  VIII,  we 
find  h = 1.77  x 10^  cm~l.  For  Qb  = 0,  the  right-hand  side  of  (47)  is 
unity,  indicating  that  c*(r,h)  = 0 in  this  limit.  This  means  that  the 
average  concentration  of  unexcited  molecules  is,  according  to  Eq.  (38) 

c = [1  - (f)3]  cT 

where  the  factor  in  square  brackets  simply  reflects  the  space  occupied 
by  the  particles,  and  can  be  taken  as  unity  for  the  range  of  a and  b of 
interest  to  us. 

The  expansion  through  terms  of  second  order  in  yQb  is  expressed 
as  a rational  function,  N/D,  with 


N = -g  (ha)  (ha)  — — (1  - ^)2(yQ b)2 

2 + — 

D = g-  (ha)  + J (f)2  + —3^  (1  + ha)  -3  J}  (yQb)2 

Further  expansion  is  not  advisable  because  the  leading  term  in  both 
expressions  is  rather  small.  The  numerator  expression  is  to  be  averaged 
according  to  Eq.  (38): 

S - [1  - (f)3]  f(ha)  +{|Ch«)(„ob)2  «f)2  - (f)3  - (f)5  + ±) 

For  our  purposes,  we  can  set  -^<<1,  £~b  at  this  stage.  Then 


N * -g(ha)  [1  + Jo  (ha)(uob)2] 
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D = ~ (ha)  + \ (1  + ha)  (y  b)2 
o J o 


In  Table  VIII  we  calculate  N/D  for  our  parameters.  For  small  particles 

or  large  R,  the  quenching  is  nearly  complete. 

To  estimate  the  effect  of  the  quenching  mechanism,  we  need  to 
* 

know  the  fraction,  f , of  excited  CO^  in  the  plume  region.  This  number 
is  not  readily  available,  and  its  estimation  involves  solution  of  non- 
equilibrium processes.  The  ratio  N/D,  calculated  above,  gives  us  the 
fraction  of  excited  CC^  molecules  that  are  de-excited  by  quenching,  based 
on  an  equilibrium  calculation.  The  excited  fraction  for  this  case  is 
(refer  to  Table  VIII) 

f eq  = /s  » \ = 4.2  x 10 

o 

which  is  a lower  limit.  In  this  limit,  it  is  clear  that  the  mechanism 
plays  a negligible  role.  The  nonequilibrium  processes  in  the  engine  com- 
bustion and  plume  formation  will  raise  f*  from  this  value.  The  change  in 
extinction  coefficient  due  to  quenching  can  be  written  as 


Aaco2  = (f)f*  CTk 


whereas,  for  h = 0,  we  will  have  the  contribution 


*co  I = (1  -f  )cTk 

<-°2  jj=0  1 


Thus,  the  fractional  change  in  extinction  coefficient  is 


a + a 
c co. 


(— ) * c k 

vD;f  T 

«c  ♦ (1  - f*)cTk 


with  ac  given  in  Eq.  (43).  The  largest  values  of  this  ratio  are  found  for 
f*  = 1,  when 
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From  Table  VIII 


16  3 

c^  = 5.0  x 10  molecules/cm 

so  c_k  =0.07.  In  Table  IX,  are  given  values  of  a and  p , for  the 
T ° c max 

parameter  ranges  of  interest.  It  should  be  emphasized  that  the  equilibrium 

"it 

value  of  f leads  to  a negligible  effect,  as  mentioned  above,  even  when  a 

c 

is  very  small,  because  the  limit  of  Eq . (57)  then  becomes 

N f*  * - 

lim  P = 7:  TT  * f N (59) 

a -»  0 01-8  D 

c 

•k 

which  is  smaller  than  f . 


Concluding  Remarks 

The  results  of  the  last  section  indicate  that  under  certain  condi- 
tions, and  in  certain  regions  of  the  (a,R)  parameter  space,  the  quenching 
effect  can  be  significant.  Th?  maximum  possible  enhancement  occurs  for  the 
largest  values  of  a and  the  smallest  values  of  R.  However,  for  this  enhance- 
ment to  be  meaningful,  the  nonequilibrium  processes  must  cause  most  of  the 
CC>2  exiting  the  jet  engine  to  be  excited.  In  the  simple  models  used  here, 
the  influence  of  nonequilibrium  processes  is  reflected  in  a reduction  in 
the  value  of  X appropriate  for  the  calculations.  The  value,  X = 500  sec  \ 
was  determined  by  using  the  measured  band  strength  of  the  v^  (4.3  pm)  line 
of  CO2,  and,  from  this,  finding  the  Einstein  "A  coefficient".  This  gives 
the  probability/second  of  spontaneous  emission  from  an  excited  C09  molecule 
in  equilibrium  with  its  radiation  field.  Collisions  are  not  explicitly 
included  in  A,  and  if  these  and  other  nonequilibrium  processes  result  in 
a higher  concentration  of  CO2*,  this  will  be  reflected  in  a smaller  effec- 
tive value  for  A (i.e.,  X). 


The  degree  of  excitation  required  in  our  model  for  useful  quenching 
is  quite  high.  The  location  of  our  typical  plume  data,  taken  from  SCORPIO  III, 
is  about  one  exhaust  diameter  back  and  one-half  of  an  exhauster  diameter  out 
radially.  At  this  point,  the  gas  has  already  cooled  to  328°K.  Unless  highly 
nonequilibrium  conditions  exist  it  does  not,  therefore,  appear  reasonable  to 
expect  a very  high  excitation  of  C^,  but  for  the  larger  values  of  a,  the 
absorption  from  carbon  is  so  reduced  that  even  the  enhancement  from  an 
only-moderately-excited  gas  will  be  a great  improvement  - perhaps  a factor 
of  50-100  in  aTQT^L*  For  the  smaller  particle  sizes,  the  quenching  enhance- 
ment is  negligible  for  all  R.  For  high  injection  rates,  only  a small  multi- 
plier may  be  introduced  by  quenching  effects,  and  except  for  the  largest  a, 
this  is  not  significant.  We  therefore  conclude  that  quenching  is  not  an 
effective  enhancement  mechanism  except  in  very  special  cases,  namely,  for 
large  a. 
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